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§ 8. Summary. 

An integral equation is obtained for the density of radiant 
energy in radiative equilibrium and solved by successive approxi¬ 
mation. The various approximations involve an arbitrary constant, 
wliich can be chosen in various ways so as to improve the 
accuracy of the solution. It is shown that the density-gradient 
taken with respect to the optical thickness ((ZB(t)/^t) has a 
logarithmic infinity at the boundary, and approximations are 
obtained for the boundary temperature T 0 in terms of the effective 
temperature T r The temperature distribution thus found is 
applied to deduce the law of darkening of a stellar disc towards 
the limb; it is shown that the value •§ for the coefficient of darken¬ 
ing is a better approximation than the Schwarzschild-Jeans value 
§; the distribution of intensity over the disc corresponds to a 
mean coefficient of about o - 6i or o'62, but there is a more sudden 
drop in intensity very close to the limb. 


Radiative Equilibrium and Spectral Distribution. 

By E. A. Milne, B.A., Fellow of Trinity College, Cambridge. 

{Communicated, by the Director, Solar Physics Observatory, Cambridge.) 

§ 1. Introductory. 

In the investigation of radiative equilibrium it is necessary to 
exercise care in assuming that the radiation will possess certain of 
the properties of black-body radiation. For example, the density 
of radiant energy for isotropic black-body radiation is aT 4 , where 
a is Stefan’s constant and T is the temperature; but it is not 
necessarily true that the energy-density in radiative equilibrium 
is equal to that of black-body radiation corresponding to the 
temperature of the matter at the goint considered. 

This point is of importance if one attempts to investigate the 
spectral distribution of the radiation.* The radiation differs from 
isotropic black-body radiation in that both the intensity and the 
spectral distribution vary from direction to direction, and the 
spectrum in any particular direction (for a given total intensity) 
will not in general be that of the corresponding black-body radia- 

* When writing this paper the author was unaware that the question of 
the spectral distribution of the light in radiative equilibrium had been 
previously discussed, but since the paper was finished his attention has been 
drawn to a paper by Lindblad, “The Distribution of Intensity in the Con¬ 
tinuous Spectra of the Sun and Fixed Stars,” Uppsala Universitets Arsskrift, 
1, 1920, which discusses very fully the theoretical spectrum for radiative 
equilibrium and the variation of colour over the disc, with detailed comparison 
with the observations of Abbot, Fowle, and Aldrich. However, it seems 
worth while making the present communication, as the method of calculation 
and the method of expressing the results are somewhat different. 
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tion. Further, the net flux, even if constant in amount, will in 
general vary in spectral distribution from place to place. Consider, 
for example, the radiation crossing a plane perpendicular to the 
direction of net flow, at a point where the temperature is T : that 
crossing from the cooler side to the hotter will he characteristic, 
on the whole, of some temperature lower than T, say T x ; that from 
the hotter side to the cooler of some temperature greater than T, 
say T 2 . If Ba(T)< 2A denotes the intensity of black-body radiation 
for temperature T between wave-lengths X and A. + dX, the spectral 
distribution of the net flow of radiation is approximately 


(B\(T 2 ) - B\(T 1 )}dX, 

or, say, 

(T, - 


which is a function of T and varies fromi point to point. The net 
flux is continually being transformed from shorter to longer wave¬ 
lengths as it advances in the direction of T decreasing. At the 
boundary of the material considered the net flux is radiated into 
space, and it is clear that its spectral distribution will not 
necessarily be that of black-body radiation of the same total 
intensity. The structure of the radiation as regards both its 
distribution in direction and its spectral distribution must he 
ascertained from detailed consideration of the laws of emission 
and absorption as expressed in the equations of transfer, with the 
appropriate boundary conditions. 

It seems worth while working out this method in detail for 
the case of the outer layers of a star. We may suppose that 
almost the whole of the energy radiated by a star is derived from 
physical processes occurring not in the thin stratum that we see 
into, but in the star’s interior; at any moment the state may he 
regarded as a steady one in which the amount of energy radiated 
externally is equal to that generated in the interior, so that the 
outer layers must take up such a state (as regards temperature- 
gradient, etc.) that they are capable of allowing the given amount 
of energy to be radiated through them. Given the total amount 
radiated per unit area, the spectral distribution of this radiation 
and the variation of this distribution with angle of emergence {i.e. 
variation of spectrum over the disc) should be determinate. 

We will first adopt the simplified treatment in which the flow 
of radiation is “linear” or “tubular,” the radiation being supposed 
confined to two streams, an inward and an outward. We will 
then modify the results by taking into account the spherical spread 
of the radiation. The first paragraph gives a few remarks on the 
general case when the absorption may depend on wave-length, but 
we shall mostly be concerned with the theoretical, spectrum for 
the case when the coefficient of absorption is the same for all 
wave-lengths. 
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§ 2. Selective Absorption and Radiative Equilibrium. 

We consider a medium stratified in planes perpendicular to the 
axis of x, extending indefinitely on the positive side of its boundary 
x = o. ' The radiation is supposed to consist of two streams parallel 
to the axis of x; let I\dX be the intensity of the outward stream 
between wave-lengths A and A + <fA at the point x, where the 
temperature is T, I \dX the inward stream. Let B\tZA be the 
corresponding intensity of black-body radiation at the temperature 
T. As in the foregoing paper (§ 3), the intensities satisfy the 
equations 

# = ^( [ »-B»). = • 0) 


whilst the equality of the energy absorbed and emitted gives 

r 00 /-oo 

2 k^B\d\= k\(I\ + I/)d\ . . . (2) 

Jo Jo ' 

Here Tc\ is the coefficient of absorption in the neighbourhood of A. 
The outward net flux of energy is 7 tFaC?A, where 


Fa = Ia-Ia' 


We shall write 

/.OO *00 /.OO *00 

I x dX = I, I I\dX = T, I B^dX — B, |Fa^A = F. 
Jo Jo Jo Jo 

From equations (1), 


( 3 ) 


whence, integrating with respect to A and using (2), 

ipl = 0, F = constant, 
dx 

as it should be. If we put 

Aa = Ia + Ia 7 ~ 2Ba, 

then 27 tAa/c, where c is the velocity of light, is the excess of 
radiant energy per unit volume of wave-length A over that of black- 
body radiation for the temperature at the point x. We can now 
write equations (1) and (2) in the form 

1 d¥\ , * , dA\ 7 T-, cfBx 

r 00 

/ k\tX\dX = o. 

Jo 

It is now clear that it is impossible for Aa to be always zero; 
for this would imply, first, that Fa = constant (from place to place), 
secondly, that &aF\oc dB\jdT, as far as variation with A is con¬ 
cerned; but since dB K /dT as a function of A alters with T and 
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therefore from place to place, these two results are contradictory, 
j Further, unless k\ is independent of A the total excess of energy 
■ A (i.e. fA\dX) is not necessarily zero. To put it loosely, 

; radiative equilibrium does not hold for the separate wave-lengths 

independently. 


§ 3. The Theoretical Spectrum for Non-selective Absorption. 

We proceed to determine the variations of I\ and 1 / with A, 
more particularly the structure of the emergent radiation I x (o), 
for the case where k\ is independent of A. Writing k for Jc\, and 
putting as usual 

f x 

| hpdx — T , 

Jo 

equations (1) become 

4 ^'-B a -Ia' ■ ■ (s) 


We saw in § 3 of the preceding paper that the corresponding 
equations for the integrated radiation of all wave-lengths have for 
their solution 


B = JF(i + 2 r) = B 0 (i + 2t) . . . (6) 


I = B + B 0 , l'-B-B 0 ... (7) 

Since B is now known, so is B x in (5); using (6) to eliminate t 
from (5), we have 

<^Ia_Ia~Bx dI/_BA —Ia' / q \ 

dB B 0 ’ dB B 0 ' W 

Solving these, 

I Ae -B/B» = _ ’ f B x er B ' B °dB + H x , 

n OJBo 

I/e-B/Bo = J. j\ e +B ' B °dB + K x ,* 

•^oJbo 


where H x , K x are independent of the temperature. Since I x ' = o 
for all A when t = o, we have K x = o. Moreover, as B->oo in the 
interior, I x /’I x '->i, for the difference between the inward and 
outward streams becomes negligible compared with the streams 
themselves. It is easily found that this gives 


Thus finally 


H x 


tt / B x e~ B/B<> (fB. 

B 0JBo 


-^o Jb 

I x ' = ^e"W B A e**dB . 
^0 Jb 0 


( 9 ) 

(10) 
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(We may also put the boundary condition at infinity in the form 
I a /B a ->i, the radiation becoming ultimately isotropic black-body 
radiation.) 

Equations (9) and (10) exhibit the spectral distribution of the 
streams of radiation at any point as derived from the temperature 
of the neighbouring material, weighted with an exponential absorp¬ 
tion factor. Since B ocT 4 , the integrals are practically integrals 
with respect to the temperature. For the stream of radiation 
escaping at the boundary we have, on setting B/B 0 = y, 

Ia(o) = ejf B A e"% . . . (n) 

Now in the usual notation 

B a oc X- 5 (efto/XET_ 1 ) — 1, 

or, putting hc/XRT 0 — a, T = T^, 

B A ocX-5(e“2/ _i -i)-i. 

Thus finally the spectral distribution of the radiation emergent at 
the boundary , in radiative equilibrium , is given by 



in which formula the only disposable constant is T 0 , which is 
itself given by the total radiation. In particular the coefficient of 
absorption k has disappeared. 

Let us compare this spectrum with that calculated on the 
assumption that the star was a perfect radiator at its own effective 
temperature T*. We should have in this case, since T 1 = 2^, 

I A (o)« 5 (e ftc A R1 ’i — 1 ) -1 /, , \ 

' ' ' 

Writing (12) in the form 

• • (it) 

we see that the two distributions could only be identical if, con¬ 
sidering the integrand in (14) as the product of a certain function 
f(y) and an exponential weighting factor e~ v , the weighted mean 
value of f(y) were /(1) for all values of a. This is in general not 
true. The function 

can be evaluated by a power series (P) in a when a < ztt, and by 
an asymptotic formula (A) when a is large.* Now the maximum 
ordinate in the spectrum is found to occur in the neighbourhood 
of a = 6, and it is such values of a that we are chiefly interested 
in ; unfortunately, for such values of a both expansions (P) and 

* See Appendix at end. 
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(A) tend to become intractable, and it is a matter of some delicacy 
to secure agreement in the overlapping values. For this reason 
values from both sources are given in the subjoined table. For 
comparison the values of 

</> 0 (a) = a. 5 (e 2 ~ ia - i) _1 

are also tabulated. The limiting values are found to be given by 
4 > i ( ol ) ^ 1 ’167a 4 , ^> 0 (a)<~ I'l89a 4 

for a very small, and by 




' a 5 e ~ 2 *“ 


for a very large. Thus 4 > 1 / 4 >o~> °'982, as a->o (large wave-lengths), 
and 0 1 /<^) o ->oo as a ->oo (small wave-lengths). 




Table I. 




i/a. 

• 

0i(a). 

A. 


<#»o(a)- 

CL» 

(P). 

(A). 

1*0 

I *oooo 

07386 


0-7585 

1 ’5 

•66 67 

2*912 


3-001 

2*0 

•5000 

7*076 


7-314 

2-5 

•4000 

I3*I24 


13-593 

3 

•3333 

2C456 


21-200 

3‘5 

•2857 

28*218 


29-317 

4 

•2500 

35'SS 1 

35-56 

36712 

4’5 

•2222 

41758 

4176 

42'920 

5 

•2000 

46-388 

46-37 

47-359 

5'5 

•l8l8 

49‘259 

49-23 

49-829 

5-8 

•1724 



50-388 

5'9°4 

•1694 



50-425 

6 

•1667 


50-36 

50-395 

6’i 

•1639 


5o'39 


6*2 

•1613 


5o-35 


6-5 

•1538 


49-91 

49-276 

7 

•I429 


48-17 

46 "8o8 

8 

•I250 


42-04 

39-300 

9 

•mi 


34’24 

30-526 

10 

•1000 


26-45 

22-291 

11 

•0909 


19-62 

• 15-483 

12 

•0833 


14-08 

10-317 

14 

•0714 


6-67 

4-148 

16 

'0625 


3-00 

1-505 

20 

•0500 


0*517 

0-1589 
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The curves of cf> t and <f > 0 are shown in fig. 1 plotted against 
1 Ja (since a oc i /A) to give a “ normal ” spectrum. It appears that, 
in spite of thdir different analytical expressions, the functions 
closely resemble one another, the curve for radiative equilibrium 
being approximately a copy of the Planck curve displaced to the 


- Theoretical spectral distribution for a star in radiative equilibrium. 

— • —• Spectral distribution for black body giving same total radiation. (Planck,) 
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violet. The maximum ordinate of the Planck curve occurs, as is 
well known, for 

2- J a = 4-965, or a= 5-904, 

whilst detailed interpolation amongst the values in Table I. shows 
that the maximum ordinate of the new curve occurs for 

a = 6-09, 

the ratio of the two being 1-031. Noticing that a oc i/AT 0 oc i/AT^ 
and recalling that Wien’s law is expressed in the form 

W T = constant, 
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we may state the difference between the curves in the form : the 
' form of the normal spectrum in radiative equilibrium is such that 
the constant occurring in Wien’s law (A,,,,.* T] = constant) must be 
; taken 3 - i per cent, smaller than for a perfect radiator, T x being 
the effective temperature of the star as deduced from the total 
radiation. In other words, we mag expect the apparent tempera¬ 
ture deduced from the total radiation to be 3'i per cent, smaller 
than that deduced from an application of the usual form of Wien’s 
law. Since in our treatment we have assumed a “ linear ” or 
“ tubular” flow of radiation, equivalent to averaging over all 
directions, we may expect the result just obtained to apply not to 
any particular portion of a stellar disc, but to the integrated 
radiation from the whole disc. 

In the case of the sun it is well known that the temperature 
deduced from the wave-length of maximum intensity exceeds that 
deduced from the solar constant, though this is often vaguely 
attributed to the “ imperfect radiating power of the sun.” The 
temperature deduced by Abbot * from the wave-length of maximum 
energy for the whole disc is 6230°, whilst that deduced from the 
solar constant is 5860°; this is about 6 per cent, lower, so that 
our correction accounts for about half of the discrepancy. But 
it is undesirable to stress this comparison with observation, owing 
to the notable departure of the solar spectrum from a black-body 
spectrum in other respects. 


§ 4. Variation of Colour over the Disc. 

The equation for the oblique radiation of wave-length A is 

cos 4 ^ = I*-Ba • • • (15) 

dr 

6 being the inclination to the outward normal, and the solution of 
this for the outward and inward streams, with the proper boundary 
conditions, is 

1\ = e TBec0 j B*(t) sec 6 e~ TBeoe dT . . ( 16 ) 

I/ = e -TSec ^| B a (t) sec ij/ e T sec ^tfr . . (17) 

J o' 

where we have put tt - 6 - iff, Ia(t, 6 ) = Ia'( t j <A) when |x<0<7r. 
Changing the variable of integration in (16) to rsec0 and putting 
r — o, we have for the spectral distribution of the radiation escap¬ 
ing at the boundary 

Ido, 0 ) = (b x (t cos 9 )e- T dr . . . (18) 

Jo 

* Abbot, The Sun { 1912), p. no; Annals Astrophys. Obs. Smithson. 
Instit., iii. (1913), P- 20 °- 
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( l 9) 


where <x= hc/XRT 0 as before. Now the total amount of radiation 
emergent in direction 6 is 

/•'JO 

(20} 


/•'JO 

I(o, 0) = f B(t cos 9)e~ T d.T 
Jo 


which would itself correspond to an effective black-body tempera¬ 
ture T\(0) given by 

/T^y 1(0, 6 ) . 

VT 0 J B 0 ' 

and hence the spectral distribution of the radiation corresponding 
to its effective temperature would be 


Tx(o, 6 ) oc 


A -5 

e a(I(o, 0)/B o )-i 


I 


(2I> 


Thus the spectrum of the radiation is not the same as the black- 
body spectrum corresponding to its effective temperature; the 
difference between (19) and (21) arises from the fact that the 
exponentially weighted mean value of B(tcos 0 ) given by (20) 
(required in (21)) is not the same as the mean value of B(t cos 0 ) 
given by (19). 

To evaluate the spectrum further, we require the explicit form 
of the temperature distribution B(t). It was shown in § 5 of the 
preceding paper that a fairly satisfactory first approximation is 
given by 

b (t) = P(£ + 2t) = B 0 (i +|t). . . (22) 


and that this gives a coefficient of darkening | 
ment with the results of further approximations, 
we have from (19) 

1 


in good agree- 
Adopting this, 


( 2 3 > 


whilst the effective temperature of the radiation would correspond 
to a black-body spectrum (from (20) and (21)) 


Ix(o, 0 )« 


A -5 


ga(i+f COS0)"£ 


I 


(24) • 


The two would only agree if the mean value of t in (23) were 
unity. 

It is now clear that the investigation of the radiation from the 
whole disc, given in § 3, really gives the case cos 6 = | (sin 6 = 745). 
Using the asymptotic formula for (23) for values of a in the neigh¬ 
bourhood of 6‘5, it was found by constructing a portion of a table 
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similar to Table I. that for cosd=i the position of maximum 
energy occurs for a = 6-51; that for (24) occurs for 

«(i +f)“ i = 4'9 6 5> or a = 6 ' 243 , 

and the ratio of these is i'043. Thus the radiation from the 
centre of the disc has*its wave-length of maximum intensity shifted 
4-3 per cent, to the violet as compared with the hlack-body 
spectrum of the same effective temperature. 

For other values of 6 , to avoid the labour of computing (23) an 
interpolation method was used. If we write 

e~ x dx 1 , . 

) 0 e a ( I +P x ~ i — 1 ~ - 1 ' ' ^ 

then £->1 asp->o, for all values of a; a; is a function of a, but 
fig. 1 suggests that it changes very slowly with a. Assuming it to 
be constant, our previous work shows that for p= 1, 

6-09(1 + ce) -i =: 4-965, or x=vz 6 , 

and for p = f-, 

6‘Si( I +§£‘)~ i = 4'96 5, or £-=1-31. 

Interpolating between the values for p — o, 1, and using the 
results in (25), we have Table II. The column headed “percent¬ 
age ” gives the percentage displacement of the position of maximum 
intensity to the violet as compared with the black-body spectrum 


of the same effective temperature. 

Table II. 

Cos 0, Sin 0 . 

Percentage. 

1 *o 

0 

4'3 

°'9 

0-436 

4 ' 1 

o*8 

o*6oo 

37 

06 

o*8oo 

2-8 

o *4 

0*915 

17 

0*2 

°*979 

o-6 

CTO 

1*000 

O’O 


The mean value of these percentages, averaged over the surface 
of the disc and weighted according to the relative brightness, is 
3*i, confirming the result obtained in § 3. 

Our final conclusion is, then, that for a star in radiative equili¬ 
brium the wave-length of maximum energy for any point of the disc 
should be displaced to the violet , relatively to the position indicated 
by the effective temperature at that point, by an amount which varies 
from zero at the limb to about 4 per cent, at the centre , the average 
for the whole disc being about 3 per cent. This shift to the violet 
as we go from limb to centre is in addition, it should be 
remembered, to that arising from the relative brightness of limb 
and centre. 
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§ 5. Comparison with Observation. 

Table III., the data of which are taken from a similar table by 
Abbot, gives the variation of colour and temperature over the 
sun’s disc ; column I. gives the temperature deduced from the 
observed wave-length of maximum intensity, using Wien’s law ; 
column II. gives the effective temperature corresponding to the 
total intensity, using Stefan’s law and Abbot’s value of the solar 
constant; column III. gives the percentage excess of I. over II., 
and column IV. the calculated excess taken from Table II. 

Table III. 


Sin e. 

Cos 0. 

^max. 

I. 

II. 

ill. 

IV. 

O '00 

1 *ooo 

0-458^ 

0 

64OO 

0 

6l60 

3‘9 

4*3 

o'55 

o ’835 

0-471 

6220 

6030 

3 ' 1 

3*8 

0-825 

o '565 

0-483 

6070 

5770 

5' 2 

2*6 

o '95 

0-312 

°‘ 5°5 

5800 

5460 

6*2 

I *2 


There is some agreement near the centre of the disc, but it is 
apparent that on the whole observation and theory are discordant. 
The percentages in columns III. and IV. change in opposite 
directions, the difference (III.-IV.) starting negative but changing 
sign and on the whole steadily increasing. The temperatures 
deduced from the wave-lengths do not decrease towards the limb 
so fast as is indicated by the total radiation. It is perhaps risky 
to draw conclusions merely from the position of maximum intensity, 
seeing that the sun’s continuous spectrum is so unlike that of 
a black body, but study of the change of distribution, towards the 
limb, of the whole spectrum confirms the result. The centre of 
the disc appears to be slightly redder than would be expected, 
and the limb bluer; the centre-limb contrast in the shorter 
wave-lengths is too little. 

This might be accounted * for by ignoring the intensity of the 
radiation and assuming a much higher value of T 0 ; this decreases 
a and consequently decreases the centre-limb contrast for any 
given wave-length (relatively to a higher value of T 0 the bluer 
wave-lengths occupy the r&le of the redder ones and the calculated 
contrast is therefore smaller). It might also be accounted for by 
consideration of the actual form of the solar spectrum, taking into 
account the varying selective absorption. Lastly, it might be 
explained as a result of scattering. Since scattering is greater for 
blue light than for red, the blue light from th» centre of the disc 
will come effectively from layers nearer the surface and therefore 
cooler than the red light, and consequently the centre of the 
disc will appear redder than it would otherwise. This has been 
pointed out by many writers ; but it does not appear to have been 
mentioned explicitly that as the centre is poorer in blue rays 

* As is done by Lindblad. 

27 
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owing to scattering, so the limb should be relatively richer in blue 
rays for the same reason, for much of the blue light scattered 
from the radial beam from the interior will tend to emerge 
tangentially near the limb. The state of affairs at the limb is 
similar to that which causes the skylight in our own atmosphere 
to be blue (especially at right angles to the sun), as Eayleigh 
showed. It must, however, be noted further that scattering 
necessitates an increased temperature gradient (for a given flux) 
to drive the net flux of energy against the increased resistance, 
and the effect of this near the centre of the disc will be partly to- 
compensate for the reddening. 


Appendix. 

When a< 27 r, the integral 

.5, 


a 


™e~ x dx 


1 e aX * “ 1 


. (26) 


may be evaluated by expanding in the power series 
a 4 J ^1 - laar* + - ^a 4 x -1 + -^a 6 x~% - . . . )x i e~ x dx (27) 


where B x , B 2 . . . are Bernoulli’s numbers, and integrating term 
by term. The coefficients involve integrals of the type 


/; 


x s e X dx, 


which are cases of the incomplete T-function; the values for s> 1 
may be made to depend on those for 1 by means of a reduction 
formula. Incomplete T-functions not being readily available, the- 
values for s = s = ^ were calculated from 

J x s e~ x dx = T(i - s) - J x~ s (^i 1 + ^1“ • • • jd x > 

the latter integral being evaluated by term-by-term integration. In 

rco 

this way the integral I x~^ s e~ s ds was calculated for s= — 1, o, 1, 3> 

. . . 39, and with these values and the values of Bernoulli’s 
numbers given by Adams,* series (27) was evaluated up to 
a=5’5; for this it was found necessary to go to 40 terms to 
obtain the result correct to four significant figures, though the 
last few terms could be estimated by extrapolation. 

For large values of a, we consider the more general integral 


a 5 e~ llP f 

Jo 


e~ x dx 


e a ( l * — 1 


0* > 


* Coll . Sci. Papers , 1, 455. 
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required in §4, which reduces to (26) when p — i. Writing 
1 4 -px — y and expanding the denominator, we have 


“ 5e "*/ 0 *i+»> -rr, ■ 09) 

which is rapidly convergent for large values of a. The terms are 
all of the type 

if e-^-^dx . . . (30) 

pJi 


for which we require an asymptotic expansion for large values of a. 
Using a well-known method, we first notice that the minimum 
value of xp - 1 + a%~ k occurs for £ = and is anc ^ 

we accordingly transform to a new variable t given by 

xp~ l + ax~ l = t 2 + 


so that t=-o corresponds to aj = (|ap)i The integral (30) now 
becomes 


I e -5p \iap)i l~ 


where 


p j dt 

FT 2 = a+y _1 - 5 P~Wap)*, 


• (30 


and the positive value of U is taken when op>4. For 
sufficiently small values of t, x may be expanded in a power series 
in the form 

* = (i a P)°( I +a 1 t + a 2 t^+ . . . ) . . (32) 

It is found that 


«! = — 


8 p 


lap 


a, 


— 3 , 


a. 


— _~-a 


640 


1 > 


"g“i j 


a 5 = 


a 9 = 


itf a i > 


3 


39 


x 64 s 


a , 


Most of the integral (31) is contributed by values of t in the 
neighbourhood of zero, and in the usual manner an asymptotic 
expansion in powers of (\ap) ~ ? is obtained by inserting expansion 
(32) in (31) and integrating term by term. The result is 


PJi 


e -xp~ l -aX 


, - 8/a, 


5 a i 


+*(»))(■+ 7f 8 “i ! ) 

+ Ke-«*( i - T VNa 1 -iL(.+N a ))] ■ ■ (33) 


as far as the term in a 4 , being as above. Here ‘h(U) denotes 
the error function 

2 

$(N) = -t- e-W. 
v 7 T J § 
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The error in stopping at any particular power of a should be 
less than the last term retained. It was desirable to check 
expansions (27) and (29) against one another for such values of 
a as were possible (for p— 1), though for values of a so small as 
4 and 5, which were practically the only ones for which a com¬ 
parison was possible, large errors would be expected in formula (33). 
It was found, however, that by stopping at the term in a"' (term 
in a 4 ), even when the next term was numerically smaller, an un¬ 
expectedly good agreement between (27) and (29) was obtained; 
the figures have been given in Table I. From a = 4toa=5'5 the 
error is never more than three units in the fourth figure, whilst if 
the next term had been retained the error would have been nearly 
a unit in the second figure. I am unable to account for this 
apparently fortuitous accuracy. In consequence of it expansion 
(33) was stopped at the term in a 4 for all the values of a con¬ 
sidered (up to a = 20). The convergence of (29) is so rapid that 
two terms sufficed when a>5. 

It is a pleasant duty to thank Professor Newall for his continued 
interest in both the foregoing investigations, and for his care and 
suggestions in revision. 

Solar Physics Observatory, Cambridge : 

1921 March i. 


Preliminary Values of Variation of Latitude at Greenwich 

during 1920. 

(Conwmnicated by the Astronomer Royal.) 

In the middle of the year 1919 the observing programme 
with the floating zenith telescope was modified, the faintest and 
brightest stars being removed from the list and other pairs inserted, 
so as to reduce the mean zenith distance of each group. From the 
observations during the period i9i9'5 to 192 ro preliminary 
places for the new pairs have been obtained and revised group 
corrections deduced. The new system of declinations has also 
been reduced to the former system. 

The period in question has unfortunately been broken by 
several long spells of cloudy weather, so that the observations are 
not distributed with the uniformity that is desirable for the 
accurate determination of group corrections. Some of the cor¬ 
rections are consequently weakly determined at present. 

With the extra material derived from a further year’s observa¬ 
tions, it is possible to improve the values given for the latitude 
variation for 1919 ( M.N ., R.A.S., 80 , 454, 1920). The following 
table gives the values of the latitude variation, in the sense mean 
minus observed latitude, during the years 1919-1920 :—* 

* In M.N ., 80 454, for “observed minus mean latitude” read “mean 
minus observed latitude. ” 
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